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We present a theory of the Raman scattering from the orbital wave excitations in manganese 
oxides. Two excitation processes of the Raman scattering are proposed. The Raman scattering 
cross section is formulated by using the pseudospin operator for orbital degree of freedom in a Mn 
ion. The Raman spectra from the orbital wave excitations are calculated and their implications in 
the recent experimental results reported in LaMnOa are discussed. 
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I. INTRODUCTION 

Since the discovery of the colossal magnetoresistance 
(CMR), much attention has been|-attracted to manganese 
oxides with perovskite structure. tTEl A variety of anoma- 
lous phenomena including gigantic decrease of the resis- 
tivity are observed in the vicinity of the phase transition 
from the charge and orbital ordered state to the forro- 
magnetic metallic one in the oxides. One of the key fac- 
tors to bring about the phenomena is the orbital degree 
of freedom in a Mn ion.a Due to the strong Hund cou- 
pling and crystalline field, two eg orbitals are degenerate 
and one of the 3d3z-2_j.2 and 3dx2-y2 orbitals is occupied 
by an electron in a Mn^+ ion. 

Extensive studies of the orbital ordering have been 
done. It is well known that orbital ordering associated 
with the Jahn- Teller (JT) type lattice distortion plays 
an important role to stabilize the layered (A) antiferro- 
magnetic ordering in LaMnOa, where spins align pacallel 
(antiparallel) in the xy plane (along the z axis).El u In 
the doped manganites, a variety of magnetic structures 
associated with the orbital orderings are reported. On 
the other hand, the dymanics of the orbital degree of 
freedom still remains to be clarified. In the orbital or- 
dered state, the collective excitation of the orbital degree 
of freedortL-tepned orbital wave, has been theoretically 
predicted. EI'Ej~tj The orbital wave corresponds to the col- 
lective electron excitation from occupied orbital to unoc- 
cupid one, i. e., the modulation of the shape of electronic 
cloud. When the orbital wave is excited, various low- 
energy properties, such^ spin wave dispersiorO, phonon 
dispersionE£l, transportEJ and thermodynamic properties, 
will be largely affected. However, the orbital wave exci- 
tations have not been observed because the experimental 
technique was limited. To clarify the dynamics of the 
orbital degree of freedom and its effects on the physi- 
cal properties in manganese oxides, it is indispensable to 
establish a method to observe the orbital wave. 

As a probe to observe the orbital wave, Ishihara et al. 
proposed the resonant inelastic x-ray scattering by whida. 
dispersion relation of the orbital wave can be detected.EHl 
One of the other possible method ta detect the orbital 
excitations is the Raman scattering.^ Although only the 



excitation at fc = and the density of states (DOS) may 
be detected, the Raman scattering has the following ad- 
vantages: its energy resolution is of the order of 1 cm~^ 
which is higher than that of the x-ray scattering and the 
different modes of the excitations are distinguished by 
the polarization analyses of the incident and scattered 
photons. Recently, the Raman scattering experiment in 
the detwinned skigle crystal of LaMnOs was carried out 
by Saitoh et alE2 They found peak structures around 
150 meV. By the detailed analyses of the polarization 
dependence and the Raman-shift energy of the spectra, 
the newly found spectra can be attributed to neither the 
multi-phonon nor the magnon excitations. Therefore, a 
theory of the Raman scattering from the remaining de- 
gree of freedom, i. e., the orbital excitation, is required 
to be developed. 

In this paper, we present a detailed theoretical frame- 
work of the Raman scattering from the orbital-wave ex- 
citations in orbital ordered manganites. We propose two 
excitation processes. One of them is analogous to that in 
the two-magnon Raman scattering in antiferromagnets. 
However, the intensities of the Raman spectra from one- 
and two-orbital wave excitations are of the same order 
of magnitudes in this process, in contrast to the magnon 
Raman scattering. In another scattering process, pho- 
ton induces exchange of electrons between Mn eg and 
O 2p orbitals, and one-orbital wave excitation is brought 
about. We formulate the scattering cross section by using 
the pseudospin operator for the orbital degree of freedom 
in Mn ions. Raman spectra from the orbital wave exci- 
tations are calculated as functions of energy in several 
polarization configurations. It is shown that the recent 
Raman experiments in LaMnOs are well explained by the 
orbital wave excitations. 

In Sec. II, the model Hamiltonian is introduced and 
dispersion relation of the orbital wave is investigated. In 
Sec. HI and Sec. IV, the excitation processes for the Ra- 
man scattering are proposed and the cross section are 
formulated, respectively. Numerical results are presented 
in Sec. V, where comparison between theory and exper- 
iment is shown. The last section is devoted to the sum- 
mary and discussion. 



1 



II. ORBITAL WAVE 

Let us set up the model Hamiltonian describing the 
electronic state of manganese oxides. We consider the 
tight-binding Hamiltonian in the cubic lattice consisting 
of Mn ions. At each site, two eg orbitals are introduced 
and t2g electrons are treated as a localized spin {St^g) 
with =3/2. We introduce three kinds of Coulomb 
interaction, i. e., the intra- (C/) and inter-(f7') orbital 
Coulomb interactions and the exchange interaction (/), 
between eg electrons at the same site. The energy split- 
ting between two eg orbitals due to the JT distortion 
of a MnOe octahedron is represented by gjT^Qx + Ql^ 
where gjx and Qi are the electron-lattice coupling con- 
stant and the normal mode of the lattice distortion, re- 
spectively. The Hund coupling (Jh) between eg and t2g 
spins and the antiferromagnetic superexchange (SE) in- 
teraction (Jaf) between nearest neighboring (NN) t2g 
spins are introduced. Among these energy parametes*. 
the intra-site Coulomb interactions are the largest .E^ 
Thus, by excluding the doubly occupied states in the eg 
orbitals, the following effective Hamiltonian for the low 
energy electronic state is derived:0 

TL — T-Lj + Hh +'Haf + 'HjT- (1) 



The first term represents the SE interaction between NN 
eg electrons given by 




where Ji and J2 are the parameters for the SE inter- 
action. Their definitions are given in Ref. 5'^ is 
the spin operator of the eg electron with S = 1/2. 
T- = cos(^mi)Tiz - sin(^m;)r„ with {m^^my^m^) = 
(1,-1,0). / denotes a direction of a bond between 
site i and site j. Ti is the pseudospin operator 
for the orbital degree of freedom defined by Ti = 

{^l'^)H-y^'a'^\^a{^)ii'dt-i'a whcrc di-^a is the annihila- 
tion operator of the eg electron at site i with spin a and 
orbital 7. This operator excludes the doubly occupied 
states of electrons. (T^^) — -f (— )l/2 corresponds to the 
state where the d^z'^_j.2 {dx2_y2) orbital is occupied by 
an electron. The second and third terms in Eq. (|l|) are 
given by 

Hh + Ti-AF = — Jh ^ Si ■ St^^i -f Jaf ^ St^^i ■ St^^j. (3) 

i (ij) 

Here, the anisotropy of the SE interaction originating 
from the tetragonal (I?4h) lattice distortion is taken into 

account as \ Ji(2)l '^1(2) ^ V '^afI "^af ~ ^0/^0 ~ 




FIG. 1. (a) The unit cell of the orthorhombic structure 
(straight lines). Open and filled circles represent Mn and O 
ions, respectively. The unit cell of the cubic perovskite struc- 
ture is shown by broken lines. Arrows indicate the displace- 
ments of O ions in the orthorhombic structures where 5i is an 
amplitude of the displacement, (b) The alternate rotations of 
MnOe octahedra around x' axis. f3 is an angle of the rotation. 

Ji(2) ~ "^1(2) ^^'^ "^AF ~ ■^AF- ^0*^^' the transfer inten- 
sity between NN d^^^ -^'^ {d^x^ -r^) orbitals along the z{x) 
axis. The last term in Eq. (|l|) is given by 

■HjT^-gjT J2 Q-i^-i- (4) 

i l—z^x 

Qii is defined by {Qiz,Qix) = Q(cos6'f^, sin6'/'^), where 
Of^ represents the mixing of two normal modes of the lat- 
tice distortion. The lattice degree of freedom is assumed 
to be frozen in this paper, since the electronic process 
which we are interested in here has much larger energy 
than the lattice excitations. 

We investigate the orbital states at the paramag- 
netic and A-A-Ej phases by applying the mean field 
approximation.^ A unit cell which includes four Mn 
sites is adopted. These Mn sites are termed Ai,A2,Bi 
and B2 (see Fig. Q). A{B) and 1(2) classify the or- 
bital and spin sublattices, respectively. As spin order 
parameters, we introduce (Siz) and {St2giz)i= 3(5'iz)), 
where (• • •) represents the thermal average. For the 
JT distortion of MnOg octahedra, C-type ordering with 
(^aT(2)'^bT(2)) = (27r/3, -27r/3) is adopted by consid- 
ering the observed lattice distortion in LaMnOs. For 
the orbital degree of freedom, we introduce the rotat- 
ing frame and adopt the order parameter as (Tij) = 
cos 9l{Tiz) + sin 9l{Tix). describes the orbital state 

at site i as |6'*) = cos ^l^-^^ — r^) + sin ^|a;^ — y"^). By 
minimizing the energy with respect to {0*}, we obtain 
(^Ai(2)' ^Bi(2)) ~ {(^Aj—dA) in both the paramagnetic 
and A-AF phases. 

The collective orbital excitations in the orbital ordered 
state are studied by applying the Holstein-Primakoff 
transformation to the pseudospin operators. Here, spins 
are assumed to be frozen. The orbital pseudospin is rep- 
resented by using the bosonic operator, ai at site i, as 
= 1/2 - ala, and f,, = 1/2(4 + In Figs. | (a) 
and (b), the dispersion relation of the orbital waves in 
the paramagnetic and A-AF states are shown. Parame- 
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ter values are chosen to be J2/J1 = 0.35, i? = 1.07 and 
9jtQ/Ji = 0.7. Ji is estimated to be about 50 meV 
from the dispersion relation of the spin wave, the Neel 
temperature for A-KF and the orbital ordering tempera- 
ture. It is shown that the orbital excitation has a gap in 
both cases. However, we note the origin of these gaps are 
different as follows: The gaps in paramagnetic and A-AF 

phases are given by 



^5jtQ(|Ji 



5J2 



VI 
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9jtQ) 



and ^ (3Ji + ^gjrQ){Ji + J2 + ^9jtQ). 

respectively, where 6a is chosen to be 7r/2. While the 
gap in the paramagnetic phase decreases with decreasing 
9JtQ, the gap in the A-AF phase remains finite. This is 
because the latter originates from the anisotropic mag- 
netic structurcu fj, = (±±') in Fig. || denotes the mode 
of the orbital wave, where ±(±') represents the relative 
phase of the Holstein-Primakoff bosons of the four orbital 
sublattices in the xy plane (along the z axis) . Among four 
modes, ( — h) is the highest and its eigenoperator is given 
by 



a[, — cosh t 
+ sinh6'(^+)(( 



(-+) 



{a-Al k — a,Bl k + a,A2 k — 0.32 k) 



.t 

Ml -fc 



'■Bl-k 



^A2 -k 



^B2-fe) 



(5) 



coshfl^ and sinhfl^ are the coefficients of the Bo- 
goliubov transformation. This eigenoperator includes the 
following linear combinations of pseudospin operators: 

„,(-+) 



^ ^ e'*^"^'(cos0A(cosh0(-+) +sinh0(-+)) 



l{unit cell) 



+e -Lb2(i)x) 

-i(cosh6'("+) - sinhe'("+)) 



~e-^^^^y+^')TB2(i)y))- (6) 



Nm is the number of sites in the unit cell and the lattice 
constants are taken to be unity. By using the Group 
theory, it is shown that the irreducible representations 

of the modes (— +), (H — ), ( ) and (++) are identified 

as Ag, Big.B^g nwA i?2g in the D2h group, respectively. 
From the irreducible representations, the allowed-mode 
symmetries of the orbital waves in several polarization 
configurations are assigned. Throughout this paper, the 
polarization configuration is denoted by (C, 77) where C 
and T] {— x,y, z,x' ,y') represent the polarization of the 
incident and scattered photons, respectively. Here, x, y 
and z axes are the directions of the bonds connecting 
NN Mn sites and x' {y') is defined by cc' = x + y {y' = 
—X + y). x', z and —y' axes correspond to a, b and c axes 
in Pnma structure, respectively, (see Fig. |^) Polarization 
configurations and the allowed-mode symmetries of the 
orbital waves are obtained as follows: 
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FIG. 2. The dispersion relation of orbital wave in (a) para- 
magnetic and (b) A-AF phases. The Brillouin zone for the 
tetragonal lattice is adopted. Parameter values are chosen to 
be J2/J1 = 0.35, i? = im,gjTQ/Ji = 0.7, and 6^'^ = 2-k/3. 
Orbital state for the paramagnetic (A-AF) phase is denoted 
by (6*^1(2), ffli{2)) = {9a, -9a) with Oa = 0.527r(0.477r). An 
asterisk in (b) shows the mode which becomes Raman ac- 
tive when the monoclinic lattice distortion is introduced (See 
Fig. @). 
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III. SCATTERING PROCESSES 

In this section, we propose the excitation processes of 
the orbital wave in the Raman scattering. Here, we con- 
sider the Raman scattering experiment carried out by 
using the visible light. It has been reported that the 
electronic energy band in the region of 0^3 eV below 
the Ferrui. level in LaMnOa consists of Mn Cg and O 2p 
orbitals.ES Therefore, it is expected that both the Mn Eg 
and O 2p orbitals play important roles in the excitation 
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processes. Taking into account these facts, the following 
two excitation processes are considered. 

One of the processes is analogous to that in the two- 
magnon Raman scattering in the antiferromagnet. Mn eg 
orbitals are considered and O 2p orbitals are integrated 
out. The schematic picture of this process is presented in 
Fig. H (a). An electron at a Mn Cg orbital is excited to one 
of NN Cg orbitals through the interaction with an incident 
photon whose energy, momentum and polarization are 
denoted as huJi, ki and Xi, respectively. By emitting a 
photon with hujf, kf and A / , one of the two electrons at 
a doubly-occupied Mn site returns to the empty Mn site. 
When the orbital states at one of the two or both the sites 
in the final state are different from those in the initial 
state, one- or two-orbital wave excitations are brought 
about, respectively. Hereafter, this process is termed the 
d-d process. It is stressed that the scattering intensities 
from one- and two-orbital wave excitations are of the 
same order of magnitudes. This is due to the transfer 
intensity between the different orbitals at NN Mn sites. 
This characteristic is highly in contrast to. the magnon 
Raman scattering in the antiferromagnet Hj ES 

In addition to the d-d process, an orbital excitation 
is brought about through the exchange of electrons be- 
tween Mn eg and O 2p orbitals. A schematic picture of 
this process is presented in Fig. |^ (b). An electron in an 
O Pa orbital is excited to the neighboring Mn eg orbital, 
where p^ represents the O 2p orbital mixing with the Mn 
eg one through the a bond. There are two electrons in 
this Mn site and one hole in an O Pa- orbital in the inter- 
mediate state. Then, one of the two electrons in this Mn 
site returns to the O site by emitting a photon. When 
the occupied orbital in the Mn site in the final state is 
different from that in the initial state, one-orbital wave 
excitation is brought about. Hereafter, this process is 
termed the d-p process. In contrast to either the d-d 
process or the magnon Raman scattering in the antifer- 
romagnet, one-orbital wave is excited in the d-p process. 



IV. RAMAN CROSS SECTION 

In this section, we present the formulation for the cross 
section in the d-d and d-p processes. 

A. d-d process 

The cross section for the orbital wave excitations is cal- 
culated in the system where the eg electrons and photons 
are coupled. We adopt the following model Hamiltonian: 

rid-d = rtg -\- tit + Hph + «e-p/i- 

The first and the second terms represent the intra-site 
electron-electron interactions and the electronic transfer 
between NN sites defined as 




initial final 
FIG. 3. The excitation processes in the (a) d-d and (b) d-p 

processes. Wavy lines represent the incident and scattered 

photons with energy uui and Uf, respectively, i and j in (a) 

represent NN Mn sites. O in (b) represents one of the 6 

oxygen p„ orbitals surrounding a Mn ion. 

Ht^'^HH + Y.^ddl-d.^a 

i'ja 

(8) 

and 

E (*?;''Crf.7'.+H.C.), (9) 

respectively, di-f^ is the annihilation operator of eg elec- 
tron at site i with spin a and orbital 7. Its energy level 
is given by Sd- tJJ is the transfer intensity between 7 or- 
bital at site i and 7' orbital at site j and is obtained by the 
second-order perturbation with respect to the electron 
transfer between Mn eg and O 2p orbitals. The explicit 

form of tj^ is given by the Slater-Koster formulas.E3 The 
third term represents the photon system as 

^p'. = EK^(4a^£a + 0. (10) 

feA 

where hj^^^ is the annihilation operator of photon with 

energy momentum k and polarization A. Finally, 

the electron-photon interaction is described by the fourth 
term in Eq. (|^) as 

ntin = -l E ^{^-V/- (11) 

Here, A{r) and j]^ are the vector potential at r = (ri -|- 
rj)/2 and the current operator, respectively. These are 
given by 
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-ikr 



ikr 



and 



h '' 



H.c. . 



(12) 



(13) 



is the unit vector along the polarization A and — 
(ri-rj). 

Being based on the Haniiltonian Tid-d, the scattering 
matrix is calculated. The initial and final states of the 
scattering are described by Ti and the intermediate state 
is described by Hf~'^. The scattering matrix in the d-d 
process is obtained as, 



Qd~d 



And 



X I 2 Ji ( -niUj 



+2J2 \ -riiUj 



Si ■ So 



Si ■ s. 




with C ~ "Ke^jih ^Juj^TJ^). E^i^f-^ represents the en- 
ergy of the initial (final) state. Ji and J2 are defined by 
Ji = I'ddtl/iU' - I + hcuj:^) + Pddtl/iU' - I - ruvf:J and 
J2 = Ild4/(U + 2Jh + K-,) + ^ddtl/iU + 2Jh - hujf:^, 
respectively, where Idd is the distance between NN Mn 
sites. Finally, by rewriting the pseudospin operators in 
Eq. (|lj) by using the Holstein-Primakoff bosons, we ob- 
tain the cross section in the d-d process as 

jd-d 4 

1 *.0„.2^2 M 



;i(27rc2 



puu' 



)S(AE-e^o))- (15) 

n^j: is the number of the boson of mode /i, momentum k 

and energy e^j^. V„f_i{k) and W,jf_i{k) are the coefficients 
of the Bogoliubov transformation connecting the boson 
operator for the v th ion to the fi th eigenmode as 



(16) 



■f^^^iy' is given by 



p _ 



j,{^^ + S,)-J2{^^-S,)\siCt. 



(17) 



with Sp = (S*, • S'.+p-). 5*^^ and CP are defined by Sp = 
sin((?* + ^rup) and CP — cos(6'* -I- ^nip), respectively. 
AE(^ hu^ —hujj:^ ) is the Raman shift energy. The cross 
section for the two-orbital wave excitation is given by 

u^.C^N 4 



i-d 



puv' 



{exf ■ P){ex^ ■ P)[K.,V,Ak)W,^{-k) 
+MP^,{k)5,y+p{V,Ak) + W,.Ak)) 

2 

x(K>'(-fc) + W^.v(-fc))] 
X (1 + ^uk)i^ + ^u'-k¥i^E er e.j:), (18) 



with 



Ji 



- J2 ( - - 5'p 



cpc:, 



-2J: 



— Sp ) CP, 



(19) 



(14) and 



MUk) 



(A( 



s. 



X SyS^f COS kp. 



(20) 



In Eq. (p^), the anti-Stokes parts are neglected for sim- 
plicity. The interaction between two orbital waves in the 
final state is neglected. This interaction is expected to 
shift the two-orbital wave Raman spectrum to the lower 
energy region as the magnon-magnoH. interaction does in 
the two-magnon Raman scattering.EZI 



B. d-p process 

The cross section from the orbital wave excitation is 
calculated by using the system where Mn Cg and O Pa 
orbitals and the electron-photon coupling are taken into 
account. We start with the following model Hamiltonian: 



(21) 



The first and second terms describe the intra-site 
electron-electron interactions and electron hopping as 
nt^ - Ht^' + i epp\s,p,sa and 



^"""^ E (^7*<.K^-+H.c.) 

17(5 cj 



(22) 



respectively, where £p is the energy level of O Pa orbital 
and piSa is the annihilation operator of the Pa electron 
at Ti -\- IpdS. Ipd is the distance between NN Mn and O 
sites. (5's are the unit vectors along x, y and z directions. 
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tjs represents the transfer intensity between NN Mn eg 
and O Pa orbitals and is given by 



.1 _ i 1 

4- _ J- I 2 2 



(23) 



7(5 



for S = x,y, z with i-^^ = — tpd is the transfer inten- 
sity between Mn d^^2_j.2 orbital at and O Pz orbital at 
'Ti + IpdZ. The electron-photon interaction is represented 
by the fourth term, 



d—p 

e—ph 



(24) 



jijS is the current operator representing the transition 
between Mn 7 orbital at fi and O Pa orbital at fi + Ip^S 
given by 



jljS = T^pd^tjS E ^IjcrPiScr + H.C. 



(25) 



The scattering matrix is obtained by the second order 
perturbation with respect to H'^Zph and is given by 

^d-p ^ ^s(Ef - E,)Y,ie,f ■ p}{e,, ■ p) J , (26) 
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FIG. 4. The Raman scattering spectra from the or- 
bital waves with (a), (b) paramagnetic and (c), (d) A-AF 
spin structures in the d-d process. is defined by 

= LujC^N/{h{2TT'f )J'} and is assumed to be indepen- 
dent of LUi and ujf, for simplicity. A value of J2/J1 is chosen 
to be 0.35. Other parameter values are the same as those in 
Fig. |. 



V. NUMERICAL RESULTS 



where p — x,y, z and 



A. d-d process 



J — Inrlt 




with A — Er 



U +IJh - A-huor, 



It is worth to mention that S'jj ^ 



eludes the linear term of Tj because the orbital excitation 
is brought about in a MnOe octahedron. Therefore, one- 
orbital wave excitation contributes to the Raman scat- 
tering. Finally, by rewriting the pseudospin operators in 
Eq. (pfl) by the Holstein-Primakoff bosons, we obtain the 
cross section in the d-p process as follows. 



jd-p 



;i(27rc2)2 
4 



— T 



X {n^fi6{AE - e^o) + (1 + n^„)5{AE + e^o)) . (28) 



The Raman scattering spectra in the d-d process are 
given by the summation of and l2~'^- Numerical 

results for the Raman spectra in the paramagnetic and 
A-AF phases are presented in Fig. |4[ Sharp spectra in 
the region 3.0 < AE/Ji < 3.6 (3.5 < AE/Ji < 4.5) 
and broad ones in the region 2.0 < AE/Ji < 7.0 
(5.0 < AE/Ji < 9.0) in the paramagnetic (A-AF) phase 
originate from one- and two-orbital wave excitations, re- 
spectively. For the one-orbital wave excitations in the 
{x,y) and {z,x) polarizations, all modes are Raman in- 
active as shown in Fig. ^ In these configurations, polar- 
ization of the scattered photon is perpendicular to that of 
the incident photon and both the polarizations are paral- 
lel to the bonds between NN Mn sites. Therefore, in these 
configurations, matrix element Si~'^ in Eq. (^4|) vanishes 
and orbital excitations are prohibited. The polarization 
dependence of the Raman spectra from two-orbital wave 
excitation are obtained as follows: 

(x, x), (z, z), (x', x') and {x' , y') active, 
{x, y) and (z, x) inactive. 

The following relations of the relative intensity are 
shown: I{x,x) ^ I{x',x') > I{z,z) » I{x\y'). These 
relations reflect the type of orbital ordering, that is, the 
C-type with 9\t^ = 0.50 77r and 0.48l7r for paramagnetic 
and A- AY phases, respectively. Due to this type of or- 
bital ordering, SE interaction between NN Cg electrons 
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in the xy plane is much stronger than that along the z 
axis. Therefore, the intensity in the (z, z) polarization is 
smaller than those in the {x, x) and {x' , x') polarizations. 
Small intensity in the (x',y') configuration is attributed 
to the interference effect between the orbital waves with 
the different symmetry. In the spectra from two-orbital 
wave excitations, several peaks and edges are shown. As 
the two-magnon Raman spectra reflect the DOS of the 
magnon, the two-orbital wave Raman spectra reflects the 
DOS of the orbital waves; position of the each peak corre- 
sponds to the van Hove singularity of DOS of the orbital 
waves. However, the spectra are not the DOS of the 
orbital waves itself because of the k dependence of the 
matrix elements M^^, . The d-d process may be usefull to 
examine DOS of the orbital waves as well as its excitation 
with k = 0. 



B. d-p process 

In Fig. H, numerical results of the Raman spectra for 
the d-p process in the paramagnetic and A-AF phases 
are shown. We find that the relative intensity of the two 
spectra from the Ag and Big modes in the {x, x) config- 
uration in the d-p process is different from that in the 
d-d process, i. e., the spectrum of the Ag mode becomes 
larger than that of the Big mode in the d-p process. As 
shown in Eq. in the Ag mode, the orbital excitation 
at each site occurs in-phase. Therefore, interference ef- 
fect in the in-phase Ag mode increases its intensity in the 
d-p process. On the other hand, in the d-d process, or- 
bital excitation is doinated by the process where the two 
electrons in the NN Mn sites are exchanged with each 
other. Because a minus sign in the scattering matrix 
arises from the exchange of electrons, interference effects 
do not occur. This reflects the AF-type orbital ordering 
in the ground state. 



C. Effect of the lattice distortion 

In the previous subsections, the Raman scattering 
spectra have been calculated in the orthorhombic crystal 
structure. However, it is reported that the lattice struc- 
ture of LaMnOa is monoclinic when-.the oxygen partial 
pressure during synthesis is reduced. tj Actually, a sample 
used in the recent Raman_scattering experiments shows 
the monoclinic structure.^ Therefore, we examine the 
effect of the monoclinic distortion on the Raman spectra 
in order to compare the present theory with the experi- 
ments. 

We adopt a model of the crystal structure which is 
schematically shown in Fig. |[ Displacements of the O 
ions indicated by arrows are taken into account up to 
the order of 0{5i/lpd) as follows. The transfer inten- 
sity between NN Mn eg orbitals in the xy plane is given 
by tl^iS) = to(0)(l + - Therefore, 



10' 


.(a) 


Big 




10° 








10' 






' \ 










::-°io-' 






'' \ ■ 


10'' 






\ ■ 


10' 








lo- 








1 


2 3 


4 






AEI 




10' 


(b) 






10° 




B.g f 




10' 












(.v,A-) j\y \ 




-^°io-' 








^ lO"* 






S 


10' 








lO' 









10' 
10° 
10' 

10' 
10" 
10 = 

lO' 
10' 
10° 
10' 
10"' 
10"' 
10"' 





3 4 

AEI J, 



3 4 

AEI J, 



FIG. 5. The Raman scattering spectra from the orbital 
waves in (a) the paramagnetic and (b)-(d) j4-AF phases in the 
d-p process. 1"^'^ is defined by /p"'' = LujC^N/{h{2-Kf)j^ 
and is assumed to be independent of uJij, for simplicity. The 
displacement of the O ions and the rotation of MnOe octahe- 
dra are set to be 5\/lpd = 0.04 and P = j^. Other parameter 
values are the same as those in Fig. ^. 




FIG. 6. The unit cell of the monoclinic structures (straight 
lines). Open and filled circles represent Mn and O ions, re- 
spectively. The unit cell of the cubic perovskite structure is 
shown by broken lines. The arrows indicate the displacement 
of O ions in the monoclinic structure where 5i^4 is the am- 
plitude of the displacement. The alternate rotations of MnOe 
octahedra around x' axis are shown in Fig. ^ (b). 
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does not change within the order of 0{5i/lpd)- 2 
does not change, either, because are proportional to 
(^Q*')^. The transfer intensity between NN Mn eg and 
O 2p orbitals changes as tpd{S) = tpd(0)(l ± 3.5 6/lpd)- 
From these facts, the dispersion relation of the orbital 
wave and the Raman spectra in the d-d process are in- 
sensitive to the oxygen displacement but the spectra in 
the d-p process are sensitive. Therefore, in the follow- 
ing, we concentrate on the d-p process. The dispersion 
relation of the orbital waves is calculated in the unit cell 
which includes 8 Mn sites. 

In Fig. 0, the numerical results for the d-p process in 
the monoclinic lattice are shown. Here, ^-AF structure 
is adopted. In addition to the spectra which appear in 
the orthorhombic lattice, new spectra marked by M are 
found. Spectral intensity of the peaks M are small com- 
pared with others. This is because the distortion from 
the orthorhombic structure is small. The new peaks at 
3.6 Ji and 4.4 Ji originate from the fact that the mir- 
ror symmetry perpendicular to—Lhe z axis is absent in 
the monoclinic structure P2i/cHj Most strikingly, there 
appears a new peak at 3.1 Ji. The peak position corre- 
sponds to the energy of the orbital wave at the X point 
marked by an asterisk in the Fig. |^ (b). This mode be- 
comes zone center mode in the Brillouin zone for the 
monoclinic lattice and becomes Raman-active. 

Now, let us compare our theoretical results with the 
experimental ones. Recently, Saitoh et al. have re- 
ported the experimental results of the Raman scatter- 
ing in a detwinned single crystal of LaMn03.E3 In addi-i 
tion to the phonon Raman spectra below 100 meVj^'O 
sharp spectra in several polarization configurations at 
120^170 meV are observed as shown in the insets of 
Fig. 0. The Raman spectra observed in such a high en- 
ergy region are usually attributed to ths multi-phonon 
excitations in transition-metal oxides.EJ However, by 
comparing the polarization and the temperature depen- 
dence of the new spectra with those of the phonon ones, 
the possibility of the multi-phonon excitation is ruled out. 
Two-magnon excitations are also ruled out because spins 
are ferromagnetically aligned in the xy plane below T/v. 
As a result, orbital wave is considered as a candidate 
for the remaining excitation. The crystal structure of 
the sample in these experiments is monoclinic_jpecause 
of the oxygen partial pressure during synthesis.Ej In this 
experiments, the 514.5 nm line (2.4 eV) of an Ar+ laser 
was used. This energy 2.4 eV corresponds to the-charge- 
transfer excitation from O 2p to Mn eg orbitalsBj There- 
fore, these experimental results are compared with our 
theoretical ones in the d-p process in the monoclinic lat- 
tice. As shown in Fig. |^, the characteristic features of 
polarization dependence and relative intensity of the ex- 
perimental Raman spectra in the region of 120~170 meV 
are well reproduced by our theoretical ones. Thus, the 
Raman spectra observed in the region of 120^170 meV in 
LaMnOa are attributed to the orbital wave excitations. 




2 3 4 5 



AE/J^ 

FIG. 7. The Raman scattering spectra from the orbital 
wave in the A-AF phase. The d-p process in the mono- 
clinic lattice distortion is assumed. The displacement of the 
O ions and the rotation of MnOe octahedra are chosen to 
be Si/lpd = 0.03,<52//pd = O.ll.^sApd = 0.07, ^4/^ = 0.15 
and /3 = j^. The anisotropy in the transfer inten- 
sity is chosen to be 7? = 1.15. Other parameter values 
are the same as those in Fig. |2| Orbital state is given 
by {0Ai(2),dBi(2),dcU2),6Di(2)) = {9 A, -0 A, A, -0 a) with 
9a ~ 0.477r. Insets show the experimental Raman spectra 
in LaMnOa at 9 K in Ref. 21. Vertical axes of the experimen- 
tal data are arbitrary. 
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VI. SUMMARY AND DISCUSSION 

In this study, we have theoreticaUy investigated the 
Raman scattering as a probe to detect the orbital wave 
excitations in orbital ordered manganites. Two excita- 
tion processes were proposed for the Raman scattering, 
i.e., the d-d and d-p processes. The d-d process is anal- 
ogous to the two-magnon Raman scattering process in 
the antiferromagnet. However, scattering intensity from 
one- and two-orbital wave excitations are of the same or- 
der of magnitude unlike the two-magnon Raman scatter- 
ing. This is due to the transfer intensity between the NN 
different orbitals. In the d-p process, photon induces an 
exchange of electrons between Mn eg and O 2p orbitals, 
and one-orbital excitations are brought about. Because 
LaMnOa is a charge-transfer type insulator where the 
optical gap is about 1 eV, the d-p process is expected to 
dominate the Raman scattering when the visible light is 
used. It was shown that the theoretical results of the Ra- 
man spectra from the one-orbital wave excitations well 
explain the experimental spectra observed in the region 
of 120-170 meV in LaMnOg. 

As mentioned in Sec. II, the orbital wave in the A-AF 
phase has a gap originating from the anisotropic spin 
structure. Therefore, we expect that the gap is sup- 
pressed by applying a magnetic field. This change will 
be reflected on the peak positions of the orbital- wave Ra- 
man spectra and will be experimentally detected. Sim- 
ilarly, the gap-less orbital wave excitation may be ob- 
served in the ferromagnetic-insulating manganites such 
as Lao.88Sro.i2Mn03 whepe the orbital ordering is ex- 
perimentally confirmcd.E3 In Sec. V, we compared the 
theoretical results for the d-p process with the recent ex- 
periments in Ref. ^ This is because the d-p process is 
dominant in the adopted energy of the incident photon. 
By changing the incident photon energy, the d-d process 
may contribute to the scattering. It is expected that the 
spectral intensities for the d-d process are enhanced with 
increasing the energy of the incident photon and resonat- 
ing with the excitation energy from the occupied lower 
Hubbard band to the unoccupied upper Hubbard one. 
When the Raman experiment using such a high-energy 
photon is carried out, we expect that the comparison be- 
tween the present theoretical results for the d-d process 
and experiments will provide much information about the 
orbital wave. It is expected that the Raman spectra from 
the two-orbital wave excitation will be observed around 
200-350 meV in LaMnOa. 
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